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Abstract 

<6 ' 

<^ • The form of the leptonic mixing matrix emerging from experiment has, in the last few years, 

O ; 

, generated a lot of interest in the so-called tribimaximal type. This form may be naturally asso- 

D \ dated with the possibility of a discrete permutation symmetry (5*3) among the three generations. 



00 
00 



Oh! 



However, trying to implement this attractive symmetry has resulted in some problems and it seems 



^v^j . to have fallen out of favor. We suggest an approach in which the 53 holds to first approximation, 

> ; 

somewhat in the manner of the old 5'[/(3) flavor symmetry of the three flavor quark model. It is 

<^ ' shown that in the case of the neutrino sector, a presently large experimentally allowed region can be 
^ '. 

■ fairly well described in this first approximation. We briefly discuss the nature of the perturbations 
\0 ■ 

i which are the analogs of the Gell-Mann Okubo perturbations but confine our attention for the 



Qh| most part to the 53 invariant model. We postulate that the 5*3 invariant mass spectrum consists of 

non zero masses for the (r, b, t) and zero masses for the other charged fermions but approximately 
degenerate masses for the three neutrinos. The mixing matrices are assumed to be trivial for the 
^ \ charged fermions but of tribimaximal type for the neutrinos in the first approximation. It is shown 

' that this can be implemented by allowing complex entries for the mass matrix and spontaneous 

breakdown of the 53 invariance of the Lagrangian. 
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I. INTRODUCTION 



It is generally considered that a full understanding of the masses and associated mixings 
of the quarks and leptons is one of the chief unsolved problems in elementary particle physics. 
While a great deal of experimental knowledge about quark masses and mixings has been 
available for quite a long time, it is only in the last few years that very detailed information on 
lepton masses and mixings has been found from a series of remarkable experiments involving 
neutrino oscillations (See for some examples, refs . These results can 

be expected to provide valuable clues toward the solution of this mass and mixing problem. 
One such clue is the fact that the leptonic mixing matrix is now known to be somewhat 



close to w 



m 
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lat is called the "tribimaximal" form; actually a number of interesting discussions 
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18| of this form have been presented over a period of years. 



For immediate convenience this form may be read from Eq. ()16|) below and it can be seen 
that one column has three equal entries while another has two entries of equal magnitude 
with the third zero. Such a structure is natural in the context of permutation symmetry 
since it is a characteristic one which brings the basis of the defining representation of (the 
permutation group on three objects) to the basis in which it is decomposed into two and 
one dimensional irreducible pieces. Of course, S3 is natural in the context of the standard 
model of elementary particles, which contains thr ee p arallel families of fermions. It has also 

Mm- 



been often discussed in the literature 
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However, it seems that this symmetry can not be exact; otherwise there would be two 
exactly degenerate neutrinos, which is not the case. It is possible that the situation may 
be analogous to another three flavor theory- the old SU{3) model in which the three fla- 
vors turned out to represent the u, d and s quarks rather than the three families we now 
know about. An initial SU{3) flavor symmetry for the then unknown strong interaction 
Hamiltonian was assumed to be broken, according to the Cell Mann Okubo hypothesis 
as: 

H = Ho + Hl (1) 

where Hq is SU{3) invariant and iJ| breaks the symmetry , leaving the iso-spin subgroup 
invariant. Actually, another symmetry breaking term, which breaks the iso-spin subgroup 
but preserves a different SU{2) subgroup is also required. A possible analogy to the 5*3 
symmetry case might be to postulate that the standard electroweak Lagrangian density has 
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a decomposition like 

£ = £o + + (2) 

where, on the right hand side, the first term is assumed to be S3 invariant, the second term 
only preserves one 5*2 subgroup of 5*3 and the third term preserves a different 5*2 subgroup 
of 5*3. There is a presumption that the first term is the dominant one. 

In the present paper we will concentrate on the fully 5*3 invariant term, Cq of Eq.©. We 
must specify what are the first approximation fermion masses and the first approximation 
mixing matrices for which we are aiming. In the cases of the charged leptons and the up 
and down quarks, the first approximation to the masses is conventionally taken to mean 
assigning the heaviest fermion of each of the three families (r, b, t) a non-zero mass and 
setting the masses of the others to zero. Considering that the quark CKM angles are all 
small it seems reasonable to consider the charged fermion mixing matrices be just the unit 
matrices in first approximation. On the other hand we would like to take the tribimaximal 
form as the first 5*3 invariant approximation to the neutrino mixing matrix. The question 
of what is a good first approximation to the spectrum of neutrino masses is perhaps not 
so clear. We will argue that an approximately degnerate neutrino spectrum is the most 
reasonable choice. Our job will be to demonstrate that the same 6*3 symmetry can give rise 
to different patterns of masses and mixings for the electrically charged fermions and the 
neutrinos. 

In this paper we will work completely in the framework of the electroweak theory, post- 
poning any possible inclusion in a grand unified theory. We will assume only a minimal 
number of fermion fields so that the neutrinos will be initially taken to be Majorana type. 
On the other hand we won't fully specify the Higgs fields in advance. We will also show that 
the model can be modified for the case of Dirac neutrinos. 

In section II we will give the needed notation for the fermion fields and shall show how the 
permutation group, S3 may act on them. Section III first contains a discussion of the form 
of the neutrino mass matrix under the assumption of S3 invariance. It is shown how it can 
be diagonalized by a matrix which is, up to a diagonal phase matrix and a possible rotation 
in a two dimensional subspace, the same as the tribimaximal form. Then, an argument is 
given that the largest part of the solution space for the three neutrino masses derived from 
the known mass differences corresponds to almost complete degeneracy. Also in section III, 
we note that this can easily be achieved by making use of the freedom to have complex 
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constants in the mass matrix. This results in a physical Majorana phase in the mixing 
matrix when the neutrinos are introduced as Majorana particles. 

Section IV contains a discussion of the mass matrix of the charged leptons when it is 
assumed that these particles acquire mass using just the standard Higgs field. There is no 
problem with the invariance allowing for a mass spectrum in which only the r lepton 
is massive, as desired. However, there is a problem with getting a trivial charged lepton 
mixing matrix, as needed to realize the tribimaximal form for the overall lepton mixing 
matrix. To solve this problem one might imagine backing up and allowing various row and 
column permutations for both the neutrino and charged lepton mixing matrices (there are 36 
X 36 possibilities) and possible two dimensional rotations in the two dimensional degenerate 
subspaces of each in the hope that the product comes out to be of tribimaximal form. This 
hope is not realized. We give a simple proof that the tribimaximal form can not be obtained 
in this manner. 

In section V, we point out a possible alternative treatment of the charged lepton masses 
and mixings which seems able to give in the first approximation, the desired trivial mixing 
matrix and non-zero mass only for the tau lepton. Instead of requiring the right handed 
charged leptons to transform under 5*3 we consider them to be singlets under S3 so that 
the relevant discrete group becomes just S^l- To construct an ^'3 invariant Yukawa term 
with the Higgs fields now requires that we introduce three different Higgs doublets, one 
associated with each fermion family. Furthermore, the Higgs potential (which may involve 
other possible Higgs fields too) should allow for spontaneous breakdown of the Ss symmetry 
in order to obtain the desired mass spectrum and trivial mixing. As an application of the 
first approximation model with three approximately degenerate neutrinos we calculate, in 
this section, the leptonic factor for neutrinoless beta decay. It depends on the value of the 
degenerate neutrino masses and the single Majorana phase in the model. 

The up and down quark masses and mixings are briefly discussed in section VI. Since we 
want to have masses only for the b and t quarks and trivial mixing in first approximation 
we use exactly the same model for the Yukawa terms as in the case of the charged leptons, 
just mentioned. Note that only the left handed quark fields, not the right handed ones, 
are assigned to transform under S3. In section VII, we briefly discuss the form of the 
characteristic matrices which reflect the invariance under any 5*2 subgroup of 5*3. These may 
be considered as plausible perturbations to give higher approximations to this model. We 
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count the number of parameters corresponding to the sum of the 5*3 matrix and either one 
perturbation matrix or two perturbation matrices invariant under different 5*2 subgroups. 
Finally, section VIII contains a brief summary and discussion. 



II. PERMUTATION SYMMETRY 



For a general orientation it may be useful to start from the part of the electroweak theory 
which contains just the minimal set of fundamental fermions and their interactions with the 
SU{2)l X f/(l) gauge fields: 



(3) 



where a=(l,2,3) are the generation indices, are the appropriate covariant derivatives for 
each term and 



pa 

GLa 



ULa 
dla 



(4) 



The notation is conventional except that we are denoting pa as the two component neutrino 
field belonging to generation a. We will need the charged current leptonic weak interaction 
terms contained in this equation: 



(5) 



Ccc = -^^iL ^Lli^Kp + h.c. H , 

where g is the weak coupling constant and W~ is the charged intermediate vector boson 
field. The hatted fields correspond to mass eigenstates and are related to the original ones 
as: 

p = Up, cl = Wcl. (6) 
Here, U and W are 3x3 unitary matrices. Finally the leptonic mixing matrix is defined as 



K = W^U. 



(7) 



It is clear that the part of the Lagrangian given in Eq.Q (which seems to be the most 
securely established) possesses a global f/(3)^ symmetry- one U{3) symmetry for each term. 
Subgroups of this very large group are candidates for extra symmetries which might be 
imposed when adding the Higgs interaction terms and potential. Discrete symmetries have 
the feature that they do not necessarily imply the existence of new Goldstone bosons when 



spontaneous breakdown to U{1)em occurs. The discrete symmetry which is staring us in the 
face is the generational permutation symmetry 5*3 and has been discussed by many people. 
Of course, there are a number of possibilities for implementing this symmetry. The simplest 
would be to have all the 5 S's's act in the same way on all the fermions. In the following we 
will discuss other possibilities too. The specific symmetry transformations take the form: 



La — SabLh, ejia — > Sab^Rhi 

where the permutation matrices S are the 6 matrices: 
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(9) 



For our purpose here it will be sufficient to use this (defining) reducible representation of 

53. 

III. NEUTRINO MASS MATRIX 



The simplest way to give masses to the three 2-component Majorana neutrino fields is to 
add to the theory, in addition to the usual complex Higgs doublet, a complex Higgs triplet: 



h 



h+/V^ h++ 

-h+/V2 



(10) 



A vacuum expectation value for will then generate the neutrino mass terms: 

1 ^ 



(11) 



To begin with. My is an arbitrary symmetric (but not necessarily real) 3x3 matrix. Requiring 
53 invariance of the Lagrangian demands invariance of this term under the transformation 
Pa^SabPb, where S stands for any of the permutation matrices in Eq. Q. This requires: 



[5,M,] = 0, 
6 



(12) 



for all six matrices S. By explicitly evaluating the commutators we obtain the solution: 
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M„ = a 
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1 



(13) 



a and (3 are, in general complex numbers while d is usually called the "democratic" matrix. 
We remark that the result of Eq. p3|) would be the same even for the case of Dirac neutrinos, 
where the initial matrix is not required to be symmetric. 

The form of Eq. (fTT|) indicates that the correct "diagonalization" of M^, should be: 



(14) 



where My is a real, diagonal matrix and f/ is a unitary matrix. Let us perform this diag- 
onalization while making contact with the historical background of the subject. It is easy 
to verify that My may be brought to diagonal (but not necessarily real) form by a real 
orthogonal matrix, R as: 



R^{al + j3d)R 



The real matrix, R is typically chosen as: 



a 

a + 3/5 
0a 



(15) 



R 



-2/V6 1/V3 
1/V6 1/V3 1/V2 
1/V3 -1/V2 



n r2 rs 



(16) 



wherein we have explicitly shown how R is constructed out of the three eigenvectors (columns 
in R) of My. In order to bring My to real diagonal form we first define the (Majorana-type) 
phases, a and r as: 



a + 3/5 = |a + 3/5|e' 



2i(7 



a = \a\e 



2iT 



(17) 



and the diagonal matrix 



P 



'-^^ 
e-'" 
e"*' 



(18) 
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Now we note that the choice 

U = RP 

yields, as desired, real positive neutrino masses: 



(19) 



U^{al + (3d)U 



\a\ 







|a + 3/?| 
lal 



(20) 



Most previous workers on models of this type have considered Mi, to be a real matrix, 
which results in a little simplification. Let us discuss this case first. It is clear from a purely 
mathematical standpoint, that the choice of matrix R has a large amount of arbitrariness. 
First, there is nothing a priori wrong with permuting the three eigenvectors, rj in Eq. ()lfj|l 
in any way. Second, there is nothing wrong with permuting the three rows of R in any way. 
In fact, there is a larger freedom of making an arbitrary rotation in the two dimensional 
subspace of the degenerate eigenvectors ri and f^. This arbitrariness leaves f2 invariant up 
to its column placement. This structure is a refiection of the possibility of decomposing the 
three dimensional representation of 5*3 into its irreducible two and one dimensional pieces. 
To illustrate the rotational freedom, we may use, instead of i?. 



R' 



(21) 



cos^ —sin^ 
s'l r2 S3 = r'l r2 rs 10 

sinC, cos^ 

Notice that r2 is left invariant under this rotation. Furthermore, the matrix element, R[^ is 
no longer zero for general ^. 

The particular choice of R written in Eq. (|16p is called the "tribimaximal" matrix. It 
seems to agree with present experimental indications when it is identified with the lepton 
mixing matrix, K in Eq. ((7j). (Of course, the identification of R with K requires the 
assumption that W is at least approximately the unit matrix.) Especially, the i?i3 element 
is known to be small and could be zero. However this model can not be exactly correct 
because it leads to two degenerate neutrino masses, which is ruled out since both — ml 
(from solar neutrino experiments) and m| — rri^ (from atmospheric neutrino experiments) 
are non-zero and not equal to each other. Thus the S3 symmetry can only be a kind of 
first approximation to Nature. Nevertheless it may be a good first approximation and hence 
interesting to study in detail. That is the point of view we will take here. 
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In Harrison and Scott have argued that, even if one accepts the predicted two-fold 
degeneracy, the 5*3 model is still not a good first approximation because the two degenerate 
masses belong to columns 1 and 3 rather than columns 1 and 2. Masses 1 and 2 are expected 
to be closer to each other because the solar neutrino mass difference is a lot smaller than the 
atmospheric neutrino mass difference. Here we would like to point out that this conclusion 
can be altered if we allow for complex values of a and /3. In that case, we can easily achieve 
a situation in which all three neutrinos are degenerate. This is a valid first approximation 
to the neutrino spectrum for a relatively large range of neutrino masses. One may notice 
this fact by adopting a typical choice |2^ of best fit neutrino mass differences: 

A = m\-m\ = 6.9 X lO'^eV'^, 
B=\ml- mil = 2.6 x IQ-^eV^. (22) 

The uncertainty in these numbers is roughly 25 per cent. From these two numbers we may- 
with a two fold ambiguity due to the unknown sign of {ml — m^- determine masses mi 
and m2 for a given choice of m^. We denote the type I solution to be the case when m^ is 
greater than mi and m2 and the type II solution to be the case when m^ is less than mi 
and m2. Note that there is an upper limit on m^ corresponding to the cosmological bound 
from structure formation jg^: 

mi + m2 + ma < 0.7 eV. (23) 

The results are shown as TableHJ parts of which have already appeared in refs. [27[. It should 
be remarked that the number of digits given for the neutrino masses are only for comparing 
them with each other. They clearly do not reflect the uncertainties of the experimental 
determinations in Eq. (j22p . 

In the large range 0.3 eV > m^ > 0.15eV, degeneracy of the three neutrino masses is a 
good first approximation to the spectrum. It is clearly less good in the lower range 0.1 eV 
> ms > 0.001 eV. Incidentally, it is seen that only in the small range just above m^ =0.0517 
eV does one have anything like a usual hierarchy. An amusing feature is that the viability of 
5*3 symmetry as a good first order approximation selects a range of favored neutrino masses. 

Actually, at first glance it may not be apparent how Eq. (|20p allows for a degenerate 
mass spectrum. To achieve an approximately degenerate mass spectrum we impose: 

\a\ = \a + 3(3\ + e, (24) 



type mi(eV) m2(eV) m^{eV) 



I 


0. 


,2955 


0. 


,2956 


0.3 


II 


0. 


,3042 


0. 


,3043 


0.3 


I 


0. 


,2446 


0, 


,2447 


0.25 


II 


0. 


,2550 


0. 


,2551 


0.25 


I 


0. 


,1932 


0, 


,1934 


0.2 


II 


0. 


,2062 


0. 


,2064 


0.2 


I 


0. 


,1408 


0. 


,1411 


0.15 


II 


0. 


,1582 


0. 


,1584 


0.15 


I 


0. 


,0856 


0. 


,0860 


0.1 


II 


0. 


,1119 


0. 


,1123 


0.1 


I 


0. 


,0305 


0, 


,0316 


0.06 


II 


0. 


,0783 


0. 


,0787 


0.06 


I 


0. 


,0000 


0. 


,0083 


0.0517 


II 


0. 


,0643 


0. 


,0648 


0.04 


II 


0. 


,0541 


0. 


,0548 


0.02 


II 


0. 


,0506 


0. 


,0512 


0.005 


II 


0. 


,0503 


0. 


,0510 


0.001 



TABLE I: Typical solutions for {mi, 1712) as nis is lowered from about the highest value which 
is experimentally reasonable. In the type I solutions is the largest mass while in the type II 
solutions ms is the smallest mass. 

where e is a real number which may be as small as one likes. A picture of this equation with 
e = 0, where a and a + 3/5 are displayed as 2 dimensional vectors making up the equal sides 
of an isosceles triangle is given in FigO The angle, ip = 2(o" — r) between them must obey 
(to first order in e): 

.znA = 34(i + ^). (25) 



It is interesting to contrast this solution with one in which the coefficient (3 in Eq. (pl 
vanishes . Then the three neutrino masses are obviously degenerate but R in Eq. (jl5p 

can be any orthogonal matrix. However, in the present case, for any finite e no matter how 
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FIG. 1: Isosceles triangle representing Eq. (|24|) in the e = limit. 

small, R can have at most a two dimensional degenerate subspace.We will see that there is 
a physical effect which persists for small e. 

Since only experiment can tell us eventually whether neutrinos are of Majorana or Dirac 
type, it seems worthwhile to show that the same results may be obtained in case the neutrinos 
are of Dirac type. In that case we have the two fields ul = p and i^/j so the mass term becomes 
—urM^ul + h.c. The diagonalization equation which replaces Eq. lfT^ is: 

UiKUL = Mu- (26) 

If, for example, both and vr transform as triplets under the same 6*3, we will again have 
the same form of My as given in Eq. (fT!^ . The bi-diagonalizing matrices then become 

Ul = RP, Ur = RP\ (27) 

where R is given in Eq. p()|l and P is given in Eq. ()18j) . This replaces Eq. ()19|l so the neutrino 
contribution to the overall lepton mixing matrix is just the same. In this case, however, 
the extra phases, P in Ul may be transformed away. There is no Majorana phase and, of 
course, no neutrinoless double beta decay. Nevertheless, the complex phase plays a role in 
the mass eigenvalues, allowing all three neutrinos to be degenerate with non zero (3. 
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IV. CHARGED LEPTON MASS MATRIX 



The charged leptons may simply obtain their masses using the conventional complex 
doublet Higgs field, 



in the Lagrangian term, 



(28) 



\f2 

C = -—eRa{M,)ab^^U + h.c. + ---, (29) 

V 

where v = \pi < 0° > and 0° = (f + 0°) / -\/2. Here, is the pre-diagonal charged lepton 
mass matrix. This matrix may, in general, be brought to a real diagonal form Mg by the 
bi-unitary transformation. 

Me = \JM,W\ (30) 

where, 

M]Me = WMlW\ MMl = UeMjUl (31) 
This leads to the mass diagonal, hatted, fields: 

Cl = Wcl, Cr = UeCR. (32) 

So far we haven't discussed the imposition of 5*3 symmetry on Mg. Because the fields 
Pa and CLa both belong to the same SU (2) ^ doublet we should, to maintain the SU (2) 
symmetry, require e^a to transform in the same manner as Pa- However e^a is not required 
to do so. Nevertheless that is the obvious first guess so let us initially require invariance of 
the Lagrangian under both transformations in Eq. (jHl). Then we must also impose 

[^,Me] = 0, (33) 

for all six permutation matrices. As before, this has the solution. 

Me = 7I + (34) 

where d is the previously defined democratic matrix while 7 and 6 are two complex numbers. 
Me will then have have the doubly degenerate real positive eigenvalue I7I as well as the non 
degenerate eigenvalue I7 + 3(5|. We define phases a' and r' by, 

7 + 3(5 = I7 + S^le^*-^' 7=|7|e2*"'. (35) 
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If these phases are used to define a diagonal matrix, P' of the form, 



e 



-it' 







P' 



e 







(36) 







e 



-it' 



we can satisfy Eq. ()30|) with the identifications. 



W = R'P' Ue = R'P'*, 



(37) 



where R' is a real orthogonal matrix of the generalized (as discussed around Eq. (PT|) ) 
tribimaximal type. We remark that it is not necessary for the ordering of the eigenvalues in 
Me (and correspondingly in Eq. ()3fi|l to be the same as in M,^. As discussed above, we can 
also permute the rows of R' and make an arbitrary rotation in the degenerate subspace. A 
reasonable first approximation to the charged lepton mass spectrum would seem to consist of 
a single massive r~ and two massless (or small mass) others (e~ and fi^). This can easily be 
achieved by setting 7 to zero and identifying the third eigenvalue to be the non-degenerate 
one. 

However, it is not so easy to get the correct lepton mixing matrix which, using Eq.(|7j). 
takes the form. 



Note that the phase matrix P'* can be eliminated by a rephasing of the field e^, which sits 
to the left of it in the physical interaction term, Eq.(jS)). The question is whether the product 
R'^R can have anything like the tribimaximal form. Ideally, one might like to keep R as 
the tribimaximal matrix and arrange for R' to approximate the unit matrix. This could be 
achieved if the complex number 6 in Eq. fl34j) could be made zero. However we have seen that 
it is necessary for 7 to be zero in order to give a realistic first approximation to the charged 
lepton mass spectrum. Clearly, having 7 = 5 = is not viable. Then, R' can not be the unit 
matrix but must be of generalized tribimaximal type. It is not expected that R' would be the 
same as R since we have seen that a different ordering of mass eigenvalues is appropriate. In 
general there are 36 x 36 discrete possibilities for the matrix R'^R corresponding to 6 possible 
row exchanges and six possible column exchanges for each and also the possibilities of two 
independent rotations in the two degenerate subspaces. At first, it seems to be a daunting 
task to study them. Nevertheless we will demonstrate a simple way to see that there are 



K = P'*R" RP. 



(38) 
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no acceptable solutions that would give a resulting K anything like the tribimaximal form 
desired. 

For this task it is convenient to represent each of the relevant matrices as a single row, 
having elements which are column vectors as in Eqs. and (j2H) . To anchor the notation 
we introduce unit vectors: 



(39) 



Now the matrices R = [ri,f2,'r^] and R' = [^1,^2,^3] may be written as row vectors with 
indexed column vector elements according to: 





1 












ni = 





, n2 = 


1 


, fis = 

















1 



fibRb 



UbRL- 



(40) 



From the first of these equations we may write He = raR^c- Substituting this into the second 
of Eqs. fHUj) we find tf = raRacR'cf, which may be more compactly written as: 



tf = faKaf. 



(41) 



Here K = KP* is the non-phase part of the lepton mixing matrix. In other words, the 
lepton mixing matrix is simply displayed as the transformation between the columns of R 
and R'. The key point for our present purpose is that one of the tf and one of the must 
be: 

^ 1 



1 

V3 



(42) 



This is seen to be the case for any generalized tribimaximal matrix which diagonalizes permu- 
tation invariant mass matrices in Eq. (PT|) . It is also clear that an arbitrary permutation of 
rows in the diagonalizing matrix does not change this vector. Without any loss of generality 
let the invariant vector in Eq. ()42|l occur as ti and r2. Then Eq. 1)411) becomes. 



h = fiKn + tiK2i + rsKsi- 



(43) 



Since the three vectors on the right hand side are linearly independent, we conclude that 
K21 = 1 and in addition Ku = K31 = 0. Similarly, the inverse relation fb = tcK~f^ yields 
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again K21 = 1 and also K22 = K23 = 0. Thus the matrix K must have the form, 



K 



X X 

1 

X X 



(44) 



where the x's stand for elements yet to be determined. Since i^' is a real orthogonal matrix 
we can finally write: 

■"O ±c ±s 

K= 1 , (45) 

-s c 

where s and c stand for the sine and cosine of some angle. The key point is that there are 
four zero elements of the matrix and one independent parameter. Therefore it represents, up 
to some possible reflections and relabeling of axes, a rotation in a two dimensional subspace. 
This structure clearly does not depend on the choices of fa and tf, we chose to identify with 
the invariant vector of Eq. ()42|) . The tribimaximal matrix of Eq. ()16|) and the generalized 
tribimaximal mixing matrix of Eg. 1)211) can not be obtained from this too simple form. 

To sum up the work so far, we have seen that S3 symmetry based on the transformations 
properties for La and e^a as given in Eq.(jHI) has some very encouraging features for a first 
approximation to lepton properties. First, the neutrino mass spectrum is consistent with 
three approximately degenerate masses. In addition, the predicted neutrino mixing matrix 
is consistent with the tribimaximal matrix, as experiment suggests. Also, the charged lepton 
mass pattern can be made consistent with a heavy t" and zero mass fi~ and e~ . Things 
would be fine if the charged lepton mixing matrix could be made consistent with the unit 
matrix so that, as in Eq.Q, the lepton mixing matrix, K could be identified with the 
tribimaximal form (up to the possibility of Majorana phases). Unfortunately, we have 
just shown that this is not possible. We will now investigate a modification of the field 
transformation properties which looks able to solve this problem. 

For clarification we remark that the problem above would not have arisen if we had not 
assumed the restriction on the S3 invariant form of the charged lepton spectrum that only 
the r particle has mass. If we could tolerate a first approximation in which (e, /i, r) are all 
degenerate, then the choice of 5 = in Eg. ()H4|1 would work. However that does not seem 
physically reasonable. 
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V. MODIFIED APPROACH FOR CHARGED LEPTONS 



As mentioned around Eq.(jHJ, there is some freedom to modify the type of the 5*3 trans- 
formations we assume. In the last section we assumed that La and e^ja transform in the 
same way so we were identifying the S^l and S^r transformations. Now let us assume that 
eRa does not transform under S^l', i-e., it transforms as a singlet under S^l- Then, in order 
to construct an invariant charged lepton mass term we need to assume that there are three 
Higgs doublets which belong to the (reducible) defining representations of 5*31,: 



where 



<ja + 0°)/v^- The charged lepton masses may then arise from a term, 



(46) 



(47) 



cM a, 



where G^^^ is a matrix of coupling constants whose form will be restricted by the S^l sym- 
metry (The S^R symmetry is not being implemented now). The charged lepton mass matrix 
is: 

{M,U = j:vaG'^l (48) 

a 

SO that the charged lepton mass term in the Lagrangian looks conventional: 



^ = - Yl (^Rc{Me)cbeLb + h.c. 



(49) 



c,6 



As before, the 5*3^ invariance demands that the three matrices, G*^'^-' have the forms: 



q{c) =^Wl + 5W(i 



(50) 



wherein the six quantities 7^'^-' and 5'^^^ are all undetermined complex numbers. Then the 
predicted form for Me is: 



t;i7(2) + A5(2) ^27(2) + A5(2) ^3^(2) + x^{2) 
t;i7(3) + A5(3) ^27(3) + A5(3) ^3^(3) + A5(3) 



(51) 



where we have denoted the sum of the three vacuum values as A = f 1 -|- f 2 -|- ^3. 
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Our goal is to see if Eq. ()5ip can yield a mass spectrum with a single massive charged 
lepton and the unit matrix for W, the charged lepton mixing matrix defined in Eq.®. This 
can be simply achieved if first, the three arbitrary constants (5'-"^ are all set to zero together 
with the two arbitrary constants 7^^-' and 7*^^^ Second, since just the bottom row remains, 
we must have a spontaneous breakdown structure in the Higgs potential which results in 
the vanishing of the expectation values vi and V2 but not ^3. The r~ mass then becomes, 
in the first approximation limit, f37^^-'. 

The Higgs potential appropriate for this model may be rather complicated since many 
additional Higgs fields could very well be present. One characteristic term illustrating how 
a basic permutation invariant quadratic form appears is: 

y = fci [$1 {eSab + (dab) ^b-k2] [$1 {eSef + C4/) $/ - fca] * + • ■ • , (52) 

where ki,k2,e and ( are constants while dab stands for elements of the democratic matrix. 

So it seems that by using, instead of a single Higgs doublet, three Higgs doublets linked 
to the three leptonic families by S^l invariance, we can obtain a suitable first approximation 
to the lepton mass spectrum and to the leptonic mixing matrix K. Specifically, the charged 
leptons have a trivial mixing matrix and only the heaviest one is massive. On the other 
hand, in this same approximation, the three neutrino masses are approximately degenerate 
and have a mixing matrix which can be taken to be of the tribimaximal form, Eq. (|T9|) or its 
generalization to include a rotation in the two dimensional degenerate subspace. 

Associated with this model in the case where the neutrinos are formulated as Majorana 
particles, is the existence of the Majorana phase ip defined in Eq. ()25p . To illustrate the 
importance of this CP violating Majorana phase we calculate the characteristic leptonic 
factor |mee| for neutrinoless double beta decay in the present model with three approximately 
degenerate neutrinos having mass, m = The general formula is: 

\mee\ = \mi{Kn)^ + m2{Ki2f + m^{Ki^)\ (53) 

where the lepton mixing matrix elements are obtained by identifying K with U given in 
Eqs. flT^ . (fTHjl and (fTBj) . This yields with e negligible, 

m 

Imp 



-- + (54) 

It is easy to check that this result would not change if we computed the matrix elements 
Kab using the rotated form of R given in Eq.()21|) rather than the one in Eq.(fTB|). Clearly 
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|mee| obeys the inequality, 

m > \mee\ > m/3. (55) 
The present experimental bound on this quantity is 

|mee| < (0.35 - 1.30)eV, (56) 

which is not much greater than the range of m for which the approximate degeneracy of all 
three neutrinos holds. Thus this kind of model may be tested in the next few years. 

VI. QUARKS 

The up and down quark mass matrices clearly have a reasonable first approximation 
in which the heaviest of each is massive while the other two have no mass. In the same 
approximation it is reasonable to have a unit mixing matrix. Then the situation for each 
is the same as for the charged leptons, as just discussed. It thus seems natural to use the 
same setup as for the charged leptons, with three Higgs fields transforming under S^l- The 
quark mass terms in the Lagrangian would take the form, 

£=-v^5:J«,^<l>li?l?% + 2v^^n^,eE€4?^2g.+ h.c, (57) 

c,b a c,b a 

which yields the up and down quark mass matrices as: 

a a 

The predictions of 5*3 invariance are 

As in the charged lepton case the first approximation to the quark mass spectrum yields, 

= Vsri^^\ rrih = V3L^^\ (60) 

while the other quark masses remain zero. For completeness we mention that the up and 
down mass matrices in the general case are brought to diagonal forms by the bi-unitary 
transformations: 

M„ = UuMuWl Md = UdMdWl (61) 
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This leads to the diagonal (hatted) states, 

ul = WuUl, ur = UuUr, di = WddL, dR = UddR, (62) 

and the hadronic part of the charged current weak interaction, 

£^^ = ^W+ULif.CdL+ h.c. + ---. (63) 

Here C = W^Wd is the quark mixing (or CKM) matrix. In the first approximation under 
consideration C = 1. 

VII. HIGHER ORDER APPROXIMATIONS 

Next let us briefly consider the perturbations which might be employed to improve the first 
approximation fermion spectra and mixing matrices. For this purpose, we first summarize 
the general form of the 5*3 invariant matrices which appear in Eqs. (|T!?|l . (jBH), (|Hn|l and 
as: 

abb 
b a b 
b b a 

This matrix contains two complex parameters. As we have seen the non-reality plays an 
important role in obtaining the proposed first approximation to the neutrino mass spectrum. 
Of course, the parameters will be different for fermions of each electric charge. The assump- 
tion of a grand unified theory would give additional relations among them. Note also that 
the 5*3 invariant matrix and possible perturbations of it directly correspond to the the mass 
matrix for the assumed Majorana neutrinos but not for the charged leptons and quarks. In 
the latter cases, the effects of the three Higgs mesons being assumed must be folded in as 
in Eqs. (l48|l and (|58|) . An additional difference is that the matrix for the Majorana neutrino 
masses must be a symmetric one. 

If we assume that the perturbations are invariant under only an 5*2 subgroup of S3, 
there are three choices. These correspond to two element subgroups containing, besides the 
identity element, S^^'^\ S^^^^ or S^'^^\ Requiring, a general matrix, M to commute with these 



(64) 
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separately gives, respectively the 5*2 invariant forms: 



'a' 


h' 


d' 




'a" 


d' 


h"' 




' d" 


d!" 


d'"' 


h' 


a' 


d 




d" 


d' 


d" 




d" 


a!" 


b'" 


d' 


d' 


d 




b" 


d' 


a" 




d" 


h'" 


a'" 



(65) 



Each of these possible perturbation matrices contains five complex parameters. In case of 
application to the Majorana neutrino masses, there will be only four complex parameters 
since we have the corresponding symmetric structures: 



'a' 


h' 


d' 




'a" 


d' 


h"' 




' d" 


d" 


d"' 


h' 


a' 


d 




d' 


d' 


d' 




d" 


a'" 


h'" 
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d 




b" 


d' 


a" 




d" 


h'" 


a'" 



(66) 



The simplest perturbation scheme would consist of the original form, Eq. (j64|) plus one of 
the three matrices in Eg .(1651) [or Eq.()66p for the Majorana neutrino case]. Each of the three 
possible sums can be seen to have five [four] complex parameters, as one can absorb the two 
S's invariant matrix parameters in the parameters of the perturbing matrix. Presumably, 
one of these three possibilities would be the best. 

The perturbation scheme which might be the closest analog of the generalized Gell-Mann 
Okubo type mentioned in the Introduction would contain the invariant part together 
with two of the matrices in Eq. ()65|) [ Eq. fjMj) ] . It can be seen that one would get the same 
number, eight [six], of complex parameters regardless of which set of two is selected. In fact, 
regardless of which set is chosen, the resulting sum of three matrices would be the same. In 
the case of the Majorana neutrino matrix, the six complex parameters already comprise the 
most general matrix. For the charged fermion matrices, nine parameters are needed for the 
most general possibility; however, all of the perturbation matrices satisfy the relation: 



Mi2 + M23 + M31 = M32 + M21 + Mi3, 



(67) 



thereby eliminating one complex parameter. Clearly, there is no point in considering a 
perturbation made as the sum of all three matrices in Eq. ()65|) . 

The most practical second approximation would seem to involve choosing just one of the 
matrices of Eq. ()65|) [Eq.()66|)]. A promising choice of perturbation, where the So subgroup 
includes the element S'*^^^-' has b 
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VIII. SUMMARY AND DISCUSSION 



For a long time the permutation invariance, Ss of the three known generations has been 
considered a very natural assumption for understanding the non-trivial spectrum of fun- 
damental fermion masses and their associated mixings. Recently, experimental data on 
neutrino oscillations have pointed to a leptonic mixing matrix which is of the so called tribi- 
maximal form; this is in essence the transformation from the natural three dimensional basis 
of S'3 to the irreducible basis and therefore can be regarded as another motivation for this 
group. We have attempted to spell out in detail a possible way in which this attractive 5*3 
permutation invariance can be consistently implemented in the standard electroweak theory. 

We started by reviewing the fact that exact 5*3 invariance is not correct and suggesting 
that the situation might be analogous to an older one where a symmetry group was very 
plausible but the exact dynamics were unknown-the SU(3) symmetry of the three quark 
model. A similar proposal is that the 5*3 invariance holds for a dominant piece but that there 
are other (presumably relatively small) pieces which preserve different 5*2 subgroups of 5*3. 
In this paper we concentrated on looking at the first approximation of exact S3 invariance. 
We showed, by comparison with results of analysis of neutrino oscillation experiments (See 
Table I) that there is a large range of allowed neutrino mass sets where near degeneracy of 
all three neutrino masses is a very good approximation. So we regarded a near degenerate 
neutrino mass spectriim together with a tribimaximal type neutrino mixing matrix as the 
5*3 invariant form to be obtained. At the same time we considered that the 5*3 invariant 
structures for the charged leptons and quarks should contains masses for the (r, t, b) states 
and zero masses for the others. The S3 invariant mixing matrices for the charged fermions 
are all considered to be just the unit matrix (no mixing to leading order). 

It is not completely trivial to obtain a model which achieves this desired spectrum. To 
get it, we included the features of i) Majorana type CP violation, ii) S3 invariance for only 
the left fermions, iii) Inclusion of three Higgs doublets linked to the three generations and 
transforming under S3 and iv) spontaneous breakdown of S3 due to only one of those three 
doublets developing a vacuum value. The motivation is given in section IV in which it is 
shown that, very generally, one can not get a tribimaximal form for the overall lepton mixing 
matrix together with the given mass spectrum of the charged leptons in the simplest way 
one would imagine. 
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A physical prediction of the model when the neutrinos are assumed to be of Majorana 
type is obtained for the neutrinoless double beta decay factor rriee already in the S3 invariant 
limit (See Eq. (j54j) ) . This gives a bound which should be testable in the relatively near future 
if the three neutrinos belong to the energy range where they are approximately degenerate. 

One natural extension of this approach would be to look at its embedding in a grand 
unified theory. The next step would be to investigate in detail the case when a perturbation 
breaking the symmetry down to an S2 subgroup is included. As discussed in Section VII, 
this would introduce two additional complex parameters for the neutrino mass matrix. 
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